Cosmic Microwave Weak lensing data as a test for the dark universe 
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Combined analyses of WMAP 3-year and ACBAR Cosmic Microwave Anisotropics angular power 
spectra have presented evidence for gravitational lensing at > 3-cr level. This signal could provide 
a relevant test for cosmology. After evaluating and confirming the statistical significance of the 
detection in light of the new WMAP 5-year data, we constrain a new parameter Al that scales the 
lensing potential such that Al — corresponds to unlensed while Al — 1 is the expected lensed 
result. We find from WMAP5+ACBAR a 2.5-a indication for a lensing contribution larger than 
expected, with Al = 3.1^'g at 95% c.L. The result is stable under the assumption of different 
templates for an additional Sunyaev-Zel'dovich foreground component or the inclusion of an extra 
background of cosmic strings. We find negligible correlation with other cosmological parameters as, 
for example, the energy density in massive neutrinos. While unknown systematics may be present, 
dark energy or modified gravity models could be responsible for the over-smoothness of the power 
spectrum. Near future data, most notably from the Planck satellite mission, will scrutinize this 
interesting possibility. 
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I. INTRODUCTION 



Results from the last decade of Cosmic Microwave 
Background (hereafter CMB) anisotropy observations 
have lead to a revolution in the field of cosmology (see 
e.g. @,@,@,|), (§,[§). Many fundamental param- 
eters of the cosmological model have now been measured 
with high accuracy. Moreover, since the standard cos- 
mological model of structure formation, based on dark 
matter, inflation and a cosmological constant, is in rea- 
sonable agreement with the current observations, CMB 
anisotropies are now considered as a cosmological labora- 
tory where fundamental theories can be tested at scales 
and energies not achievable on earth. 

One crucial test concerns the nature of the dark energy 
component and the validity of General Relativity (GR, 
hereafter). The simple fact that supernovae type la ob- 
servations are in agreement with an accelerating universe, 
which is puzzling in several theoretical respects, calls for 
the deepest possible investigation of dark energy and for 
a continuous test of GR. 
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CMB anisotropics are mainly formed at redshift z ~ 
1000 when cither dark energy or modifications to GR 
appear to be negligible. However, while CMB photons 
travel to us, they are affected and distorted by other, low 
redshift, mechanisms, that could help in understanding 
the nature of the accelerating universe. 

The so-called late Integrated Sachs- Wolfe effect, for 
example, generated by the time- variation of the gravita- 
tional potential field along the CMB photon's line of sight 
in dark energy dominated universes, has already been 
detected by more than five groups by cross correlating 
galaxy surveys with anisotropies at very large angular 
scales (see e.g. ||). While the statistical significance 
of the effect is still under 5c, the detection represents a 
crucial test for dark energy |7j. 

On scales of ten arcminutes and smaller, the interac- 
tion of the CMB photons with the local universe starts 
to be dominant with second order anisotropies arising 
from weak lensing or scattering of the CMB photons off 
ionized gas in clusters and large scale structure (Sunyaev- 
Zel'dovich- SZ effect). 

Weak lensing of CMB anisotropics could provide use- 
ful cosmological information. Gravitational lensing can- 
not change the gross distribution of primary CMB 
anisotropies, but it may redistribute power and smooth 
the acoustic oscillations in the CMB power spectrum (sec 
e.g. Only in the tails of Silk damping (||, at 

I > 3000) the lensing contribution start to change the 
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power spectrum significantly. Higher signal-to-noise can 
be achieved by correlating power in different directions 
on the sky, effectively using the four-point function signa- 
ture imprinted by lensing to reconstruct the line-of-sight 
integrated matter distribution 1 . 

The strength of the weak lensing smoothing is related 
to the growth rate and amplitude of the dark matter 
fluctuations. Since both dark energy or modified gravity 
significantly affects these perturbations, a measurement 
of the CMB lensing, through its high-£ smoothing, can 
in principle be a useful cosmological test (see e.g. p0|). 

The recent claim made by the ACBAR collaboration 
(fill) for a detection of weak lensing, based solely on 
smoothing of the angular power spectrum, opens the 
opportunity for this kind of analysis. To first order, 
lensing causes the primordial peak structure to be less 
pronounced, as gravitational potential fluctuations on 
large scales mix the various scales in the primordial 
CMB power. Based on the effect on the power spec- 
trum, the ACBAR collaboration has reported a A% 2 = 
9.46 between the lensed and unlensed best fits to the 
WMAP+ACBAR data, which translates into a > 3cr de- 
tection of CMB lensing. 

In this paper we further analyze this result and we 
study the possible cosmological implications. In the next 
section we phenomenologically uncouple weak lensing 
from primary anisotropies by introducing a new param- 
eter Al that scales the gravitational potential in a way 
such that Al = 1 corresponds to the expected weak lens- 
ing scenario. We then constrain this parameter with cur- 
rent CMB data, we evaluate the consistency with Al = 1, 
the correlation with other parameters and with other sys- 
tematics such as SZ. We will report a ~ 2cr preference 
for values of Al > 1. We will then discuss some possi- 
ble cosmological mechanisms that can increase the CMB 
smoothing, namely an extra background of cosmic strings 
and modified gravity. 



II. ANALYSIS METHOD 

Weak lensing of the CMB anisotropies enters as a con- 
volution of the unlensed temperature spectrum Ci with 
the lensing potential power spectrum Cf (see [||). This 
convolution serves to smooth out the main peaks in the 
unlensed spectrum, which is the main qualitative effect 
on the power spectrum on scales larger than the ACBAR 
beam, or 6'. 

The weak lensing parameter is defined as a fudge scal- 
ing parameter affecting the lensing potential power spec- 
trum: 



Cf ->A L C?. 



(1) 




1 This type of estimator has recent- 
order 3 — a in the WMAP data J 
galaxy surveys. 
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FIG. 1: This figure shows the effect of varying Al parame- 
ter. The curves with increasingly smoothed peak structure 
correspond to values of Al of 0,1,3,6,9. 



In other words, parameter Al effectively multiplies 
the matter power lensing the CMB by a known factor. 
Al = is therefore equivalent to a theory that ignores 
lensing of the CMB, while Al = 1 gives the standard 
lensed theory. Since at the scales of interest the main 
effect of lensing is purely to smooth peaks in the data, 
Al can also be seen as a fudge parameter controlling the 
amount of smoothing of the peaks. The Figure illus- 
trates this effect of varying Al on a concordance cosmo- 
logical model. 

In what follows we provide constraints on Al by an- 
alyzing a large set of recent cosmological data. The 
method we adopt is based on the publicly available 
Markov Chain Monte Carlo package cosmomc flTj with 
a convergence diagnostics done through the Gelman 
and Rubin statistics. We sample the following eight- 
dimensional set of cosmological parameters, adopting flat 
priors on them: the baryon and cold dark matter den- 
sities <x>b and u) c , the ratio of the sound horizon to the 
angular diameter distance at decoupling, 8 S , the scalar 
spectral index ns, the overall normalization of the spec- 
trum A at k = 0.002 Mpc" 1 , the optical depth to reion- 
ization, r. Furthermore, we consider purely adiabatic 
initial conditions and we impose spatial flatness. We also 
consider the possibility of a massive neutrino component 
with fraction f v > and, finally, we add the weak lensing 
parameter Al- 

Our basis data set is the three-year WMAP data ||] 
(temperature and polarization) with the routine for com- 
puting the likelihood supplied by the WMAP team. As 
we were approaching complctition of this paper, the five 
year WMAP result data became available (||], ||). We 
have therefore checked that our results are stable with 
respect to the new data. 

We add the high quality and the fine-scale measure- 
ments from the ACBAR experiment (|ll[) by using the 



data set provided by the team, including normalization 
and beam uncertainties, window functions and the full 
error covariance matrix. 

Finally, we also consider an "everything" data set. 
This adds other CMB experiments Boomerang 2K2 
(g§), CBI (0), VSAE (pi), the large scale struc- 
ture data in form of Red Luminous Galaxies power spec- 
trum (| 12 1) and the supcrnovae measurements from SNLS 
(fl6|), a prior on the Hubblc's constant from the Hubble 
Key project (|1^|) and, finally, a Big Bang Nucleosynthe- 
sis prior of uj b = 0.022 ± 0.002 at 68% c.l. to help break 
degeneracies. 



III. BASIC CLAIM AND ITS STATISTICAL 
SIGNIFICANCE 

First we run two sets of Markov-chains with Al fixed 
to or 1. We measure the difference between the best fit 
lenscd model and the best fit unlensed model of A% 2 = 
9.34, which is in excellent agreement with the original 
claim by the ACBAR team (A\ 2 = 9.46). Since both 
models have the same number of degrees of freedom, this 
has been interpreted in |ll| as > 3<7 detection of the 
lensing signal. 

Can this difference be attributed to a single point? As 
can be seen in the Table H where we report the contribu- 
tion to the overall \ 2 coming from the individual points 
(using the full covariance information) the answer is neg- 
ative: the difference appears as randomly distributed 
across the 26 ACBAR points. 

The effect is also marginally present in the WMAP 
third year and five year data. Considering only the 
WMAP third year result we found a A% 2 ~ 1.6 between 
the Al = 1 and Al = maximum likelihood model. 
Considering the newly released WMAP five year data 
(® U) which extend to higher I we get A% 2 ^3.1. 

We can ask the question of significance in the Bayesian 
way, which should be more accurate in this relatively low 
signal-to-noise regime. In the Bayesian theory, the rela- 
tive probability of a model (assuming the prior probabil- 
ities on each model are the same to start with) is given 
by its evidence, which is the integral of likelihood over 
the prior (see e.g. g§, §1). 



E = J L(6)d N 9 



As shown in fl2C|], the evidence can be written as 



( j 



(2) 



(3) 



where L max is the likelihood at the most likely point and 
Vl and Vn are suitably defined volumes of posterior and 
prior. 

The crucial point for this paper is that the evidence 
ratio for the lenscd and unlensed model can be written 
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TABLE I: This is the contribution to the overall \ 2 coming 
from the individual points, using the full covariance informa- 
tion. This quantity is not constraint to be positive, as it is 
equal to A\ 2 = ((d — i) T C~ 1 )i(d — where d denotes data 
vector, t denotes theory vector and C is the covariance ma- 
trix and there is no summation over repeated indices. This 
table shows that there are no significant outliers in the data 
as the overall contribution to x 2 is evenly distributed across 
the bins. The signal is coming from a range of scales. 



simply as 



A log E = A log L„ 



AVr 



(4) 



since the prior volumes cancel exactly for the same un- 
derlying parameter space. The posterior volume can be 
roughly estimated as 



(5) 



where Ui are the marginalized estimates of the errors 
from the Markov Chains. A considerably better estimate 
would be to take the full error covariance into account, 
however, the models are so close that the noise in esti- 
mating the error covariance would probably dominate. 
This allows to estimate the evidence ratio to be 
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4.67 + 0.075 = 4.75 



(6) 
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WMAP3 + ACBAR +strings 2.9; 
WMAP3 + ACBAR +SZ1 3.1; 
WMAP3 + ACBAR +SZ2 3.0; 



3+2.3 
2-1.8 
0+2.2 
0-2.0 
0+2.3 
0-1.8 



TABLE II: . This table shows results for constraints on the 
Al parameter. We report one and two sigma errors. Note 
that all results are statistically compatible with the standard 
prediction of Al = 1 at the level of 2-3 a. 



FIG. 2: This figure shows the ACBAR data with Ce spectrum 
predictions suitably multiplied to show the structure of the 
peaks more clearly. 

The net result is that the evidence difference is dom- 
inated by the best-fit effect: both theories are equally 
good at fitting the available parameter volume, however, 
the best-fit model is considerably better for the lensed 
model. In fact, the volume factor strengthens rather than 
weakens the evidence for lensing in the ACBAR data. 

IV. VARYING A L 

However, the anticipated forecast for the ACBAR de- 
tection from Fisher matrix analysis is only at about 1- 
sigma level. How are the ACBAR results at a so much 
higher confidence limit? 

The Figure || show the ACBAR points plotted against 
C(£(£ + l)/27rexp(^/500), where the exponent has been 
chosen to roughly counter-act the Silk's damping. We 
see that there is a weak "chi-by-eye" evidence that the 
ACBAR data are actually overly-smooth given the theo- 
retical predictions and that this over-smoothness is driv- 
ing up the detection. 

We have therefore performed additional runs where 
we let A L vary. We consider the case with 
WMAP3 data alone, with WMAP3+ACBAR data, and 
WMAP3+everything data sets. Our results are summa- 
rized in the top half of Table [il] and in Figure ^. 

We see that the results prefer values of A L which are 
considerably higher than unity. As we show below, the 
result is not affected by the inclusion of the Sunyacv- 
Zeldovich component. Therefore, the detection is com- 
ing from the smoothness of peaks, rather than excess of 
power on the smallest scales. This can also be seen "by 
eye" in the Figure [|. 

The level of confidence for excess is above 2a (except 
for the WMAP data alone case which is ~ la) but less 
than three sigma away from one. In agreement with a 
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FIG. 3: Marginalized 1-D likelihood disribution for Al 
for different datasets considered: WMAP3-alone (solid 
bold), WMAP3+ ACBAR (dotted), WMAP3+" everything" 
(dotted bold), WMAP3+ACBAR+strings (solid), 
WMAP3+ ACB AR+SZ 1 (dashed), WMAP3+ACBAR+SZ2 
(dotted-dash). 

simple Fisher matrix forecast, we find a standard devi- 
ation of the lensing amplitude of AA L = 1. We also 
looked for correlations between Al and other parameters 
and found them to be negligible for all other parameters. 

Also in Table we report a similar analysis but now 
cosidering the recent WMAP 5-year data release. As we 
can see, while the error bars arc slightly reduced, the 
new data confirm the results obtained with the previous 
WMAP 3-year data. 

How shall we interpret these results? Let us consider 
three possibilities: 

1 . The result is a statistical fluctuation. We note that 
the result is less than three sigma away from the 
theoretically most expected value of 1. The sim- 
plest explanation is that this is 2-3 a statistical 
fluctuation, with data fundamentally in agreement 
with the lensed CMB theory. However, at the same 
time, the unlensed theory is deep in the tails of the 
Al probability distribution and therefore has a con- 
siderably worse x 2 - In other words, ACBAR had a 
lucky noise realization to be able to claim detection 
of lensing. 



5 



2. Hint of new physics. It is possible that new physics 
is responsible for over-smoothness of the power 
spectrum. This is obviously the most interesting 
option. We explore these possibilities in further 
detail in the following two sections. 

3. Unknown foregrounds or experimental systematics. 
A natural possibility is an unaccounted systematic 
in the experiment itself. CMB experiments are 
intrinsically difficult and despite many jack-knife 
tests that the authors have performed one should 
not exclude a possibility of a systematic that has 
slipped through. We discuss in the next section the 
possibility of an unknown foreground component. 



V. ADDITIONAL COMPONENTS 

We will now consider whether there could be an addi- 
tional component that could bring about smoothing. It is 
possible that a smooth continuous component could lead 
to an effective smearing of the peaks when the adiabatic 
component where reduced by an appropriate amount. In 
order to check this idea we have tried to add three differ- 
ent templates, whose amplitude was allowed to be free- 
floating: 

• SZ template I. A template expected from the 
Sunyacv-Zcrdovich effect as given by the analytic 
model of Komatsu and Seljak (pl|). 

• SZ template II. A similar template based on 
smoothed particle hydro-dynamics simulations p7i . 

• String template. A template corresponding to "wig- 
gly strings" of [^2|. Note that the exact shape of 
the strings corresponding to a particular model is 
unimportant. The basic question we try to address 
is if a broad, featureless addition to the power spec- 
trum can bring about a sufficient change. 

Effect of these templates on the value of Al is very 
small as shown in the results in the Table |l| and Figure ^. 
We conclude that while the data allow for some amount 
of extra smooth component, it by no means changed the 
"detection" of lensing. 



VI. NON STANDARD MODELS 

It is certainly important to investigate if there is any 
possibility to explain the anomaly through a mechanism 
based on non-standard physics. As we pointed out in 
the introduction both dark energy and modified gravity 
can change the growth and amplitude of dark matter 
perturbations and thus enhance in principle the CMB 
weak lensing signal. 

Dark energy could affect the growth by changing the 
expansion history and by gravitational feedback of the 



per turbations in the dark energy component (see e.g p3| , 
fl24j ). However quintessence scalar field models are gener- 
ally unable to produce deviations larger than few percent 
of the CMB weak lensing signal. More exotic dark energy 
models with non-zero anisotropic stresses (see e.g. p5| , 
p6| ) could be responsible for the anomaly. 

One should however consider the possibility that grav- 
ity is more complicated than anticipated by Einstein and 
that this modification causes more lensing. A feature 
common to a broad range of modified gravity theories is 
a decoupling of the perturbed Newtonian-gauge gravita- 
tional potentials <f> and ijj. Whereas GR predicts ip = <fi 
in the presence of non-relativistic matter, a gravitational 
slip, defined as ip ^ 0, gencrically occurs in modified 
gravity theories (see e.g. @, ||, || |§ ||, || Hfl). 



Gravitational lensing phenomena depend directly on 
the sum of the two gravitational potentials and is strongly 
affected by 
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l Dy a gravitational slip (see e.g. 
|38| , |39| , M). It is therefore interesting to investigate if 
Al > 1 could be explained with modified gravity and to 
more quantitatively connect this parameter to modified 
gravity theories. 

Since a very large number of models have been con- 
ceived here we use the parametrization of Daniel ct al. 
2008 (|^l|), which is simple and easy to apply to several 
models. In this parameterization the gravitational slip is 
given by a function w(z) such that i\) = (1 + zu)(f> and is 
parameterized by a single parameter tjjq defined as 



^o^(l + z)- 3 - 



(7) 



i.e. it starts to be relevant at dark energy (or modified 
gravity) appearance. 

Following jll|, we can easily approximate the relation 
between Al and zu as 



A L (w) 



\Gacdm(z = 2) 



2 + w 



(8) 



The difference in growth factors is evaluated at z = 2, 
since the lensing kernel peaks at that redshift. Larger 
values of vjq correspond to larger values of Al- A value 
of n7o ~ 1.5 could produces very similar results on the 
CMB to Al > 1.5 and thus bringing the signal inside the 
1 — a cl. According to |IlJ] this range of values of zuq is 
in agreement with the measured temperature anisotropy 
signal on very large angular scales but is at odds with 
the recent ISW detections. 



VII. SYSTEMATICS 

Let us in this section investigate what kind of system- 
atic effect could mimic the observed over-smoothing in 
the data. As we have shown in Table | the effect is not 
coming for a particular rogue data point or a small range 
of scales. This further constrains possible sources. 
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First we note that most effects that produce smooth- 
ing in real space, such as inaccurate characterization of 
the beam or pointing will induce multiplication of the 
real power spectrum by the Fourier transform of the ef- 
fective beam. This is unlikely to produce the additional 
smoothing required to explain the hint of an anomaly 2 . 

Atmospheric fluctuations could play a role. However in 
this case the effect would appear as an additional smooth 
background component and, as shown in Table ||, our 
result appears stable under this assumption. 

It may however be possible that an unaccounted for 
systematic is present in the data set provided by the 
ACBAR team, especially in the assessment of the sky 
window functions. Sky coverage of the ACBAR telescope 
is very complicated pattern of many fields with somewhat 
fuzzy edges. A poor characterization of the variation of 
noise across the fields could, in principle, lead to to the 
effect observed here. It however exceeds the scope of this 
paper to investigate this thoroughly. 

Finally, it is possible that the error has been induced 
in the final power-spectrum estimation step of the data- 
reduction procedure. The maximum-likelihood estimator 
employed by the ACBAR team in principle assumes a 
step- wise power spectrum and the real shape of the power 
spectrum has to be accounted for carefully, especially at 
the signal-to-noise present in the ACBAR data. 

It is clear that at the present stage systematic effects 
can not be ruled out and more data is needed. Fortu- 
nately, weak lensing will also produce a B-mode polariza- 
tion signal that, if observed, will provide a fundamental 
cross-check. 



VIII. CONCLUSIONS 

We have reanalyzed the ACBAR angular power spec- 
trum in light of the recent detection of a lensing signal in 
their angular power spectra. We tracked this down to a 
hint of over-smoothness in the power spectrum, detected 
at ~ 2.5a statistical significance. This over-smoothness 
pushed the theory without lensing deep into the tails and 
make it a poor fit to the data. 

If interpreted as real, there are several interesting pos- 
sibilities. A modified gravity can induce extra amount of 
lensing and we show that a gravitational slip could bring 
the discrepancy to sub 1-a level. 

How does this compare with other detection of lens- 
ing in the CMB. Two groups (E3],|43 ) have searched for 
CMB lensing by correlating WMAP data. The WMAP 



data have lower intrinsic potential for measuring CMB 
lensing than ACBAR, however by using more informa- 
tion than the smearing of the Ci structure (i.e. an op- 
timal quadratic estimator), and by correlating to galaxy 
surveys, they were able to find significant evidence at the 
3-er level. While the mean value found is close to unity, 
these previous results allow considerable freedom in over- 
all amplitude and a reasonable fit can be obtained with 
values of Al lying somewhere in between. In particular 
Al ~ 1.7 is compatible with both probes at less than 2 
standard deviations. However a possible interpretation 
is that lensing is somehow enhanced inside the ACBAR 
field of view, which is only 1% of that of WMAP. It will 
be very interesting to apply quadratic estimator tech- 
niques using the full four-point function information to 
the ACBAR maps . As the statistical error (based on 
Fisher matrix forecasting) for this probe is about 4 times 
smaller as compared to the smearing of acoustic peaks 
investigated here, we anticipate that this will shed light 
on the findings of the current paper. 

Looking at closer measurements of lensing, the weak 
lensing tends to give values of a% that seem only 
marginally higher than that of WMAP3 (see for example 
pi) . Q) and consistent with the more recent WMAP5 
measurements These measurements would limit the 
value Al < 1.2. However, the redshift spans involved are 
considerably smaller with typical redshifts probed being 
around ~ 0.5. Therefore, the drastically different source 
redshifts imply that these results are not in direct contra- 
diction and that it is conceivable that modified gravity 
models can be constructed that satisfy all observational 
constraints. 

Maybe less excitingly, but more realistically, the fea- 
ture should be interpreted as a noise realization fluctua- 
tion or explained by unaccounted systematics. 

Future experiments as Planck, especially with the help 
of polarization data, will soon shed light on this intrigu- 
ing result. 
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